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Abstract —The mixed radio frequency (RF)/free-space optical 
(FSO) relaying is a promising technology for coverage im¬ 
provement, while there lacks unified expressions to describe 
its performance. In this paper, a unified performance analysis 
framework of a dual-hop relay system over asymmetric RF/FSO 
links is presented. More specifically, we consider the RF link 
follows generalized or rj-n distributions, while the FSO link 
experiences the gamma-gamma distribution, respectively. Novel 
analytical expressions of the probability density function and 
cumulative distribution function are derived. We then capitalize 
on these results to provide new exact analytical expressions of the 
outage probability and bit error rate (BER). Furthermore, the 
outage probability for high signal-to-noise ratios and the BER 
for different modulation schemes are deduced to provide useful 
insights into the impact of system and channel parameters of 
the overall system performance. These accurate expressions are 
general, since they correspond to generalized fading in the RF 
link and account for pointing errors, atmospheric turbulence and 
different modulation schemes in the FSO link. The links between 
derived results and previous results are presented. Finally, 
numerical and Monte-Carlo simulation results are provided to 
demonstrate the validity of the proposed unified expressions. 

Index Terms —Free-space optical communications, outage 
probability, bit error rate, atmospheric turbulence. 

1. Introduction 

Due to the merits of free license, low cost and high 
bandwidth, free-space optical (FSO) communication is be¬ 
coming one of the promising technologies for indoor and 
outdoor wireless applications. Some typical applications of 
FSO systems include “last mile” access, indoor positioning, 
disaster recovery, military applications, underwater system, 
device-to-device communications, and video transmission, etc. 
ffl-Q. On the other hand, dual-hop relaying has been widely 
adopted in the context of radio frequency (RF) communication 
systems because of the extended coverage area and enhanced 
receive signal induced by spatial diversity. In traditional RF/RF 
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dual-hop communication systems, the scarcity of licensed 
spectrum in conjunction with fast growing demand of high 
data rate is one of the most crucial limitations. To this end, 
the FSO communication has been considered to be used in the 
relaying as a complement to RF counterpart due to its desirable 
features: First, the RF and FSO links operate on different 
frequency bands to avoid significant interference in RF/RF 
relay systems; Second, the FSO link can achieve maximum 
system capacity by aggregating multiple users’ information. 

Therefore, the mixed RF/FSO dual-hop relay system has 
been recently proposed a-in and refers to the case when 
RF transmission is used at one hop and FSO transmission 
at the other. Such topology is quite different from hybrid 
RF/FSO system, which uses parallel RF and FSO links for the 
same path. In IH, a performance analysis of a dual-hop relay 
system composed of mixed RF/FSO links was first conducted. 
The impact of pointing errors on such topologies has been 
analyzed in ISl. The authors in O derived an exact closed- 
form expression of the end-to-end outage probability of the 
mixed RF/FSO system. Moreover, the performance of dual¬ 
hop systems with Nakagami-m and gamma-gamma (FF) or 
double GG O fading channels was investigated in (Tl, O. 

Although the performance of mixed RF/FSO relaying sys¬ 
tems over fading channels has been extensively evaluated in 
terms of outage probability and error rate, most of existing 
pioneering studies simply assume the Rayleigh or Nakagami- 
m fading channels for the RF link. While such assumption 
extensively simplifies some mathematical manipulations, these 
distributions fall short of capturing the actual fading statistics 
of the RF link, because their tail does not seem to yield a 
good fit to experimental data HOI . Motivated by this fact, 
in this work we examine the performance of a dual-hop 
RF/FSO transmission system experiencing more general n- 
/i/FF and //-/x/FF fading conditions. The 7^-/x and n-ji fading 
were proposed in ifTll and has been verified to provide an 
excellent fit to experimental data. Moreover, it provides a 
unified framework for existing RF channel models, since it 
regards the Rayleigh (j] —> 0,/x = 1), Rician (/x = 1) 

and Nakagami-m (r] 0 ,k, 0,/x = m) distributions as 

special cases. For the FSO link, we use the FF distribution 
to describe the weak or strong atmospheric turbulence fading, 
since the FF distribution has simpler statistical characteristics 
than the double GG distribution, and has been proved to 
provide a good agreement between the theoretical result and 
the corresponding experimental data ifT^ - lfT^ . Moreover, the 


2 



Fig. 1. The mixed RF/FSO dual-hop relay system with S, R and D denoting 
the source node, relay node and destination node, respectively. Moreover, the 
first hop (RF link) experiences a K-fi or ri-fi fading, and the second hop (FSO 
link) follows a FF fading. 


pointing errors, which is caused by an unavoidable alignment 
between transmitter and receiver due to thermal expansion, 
dynamic wind loads, and weak earthquakes, is one of the major 
issues leading to a severe performance degradation of the FSO 
link (TSl, lO. Therefore, we study the effect of pointing errors 
on the performance of dual-hop RF/FSO systems. 

The contributions of this paper are summarized as follows: 

• Novel analytical expressions of the probability density 
function (PDF) and the cumulative distribution function 
(CDF) for the mixed RF/FSO dual-hop relaying system 
are derived in terms of Meijer’s G-function, which can 
be easily evaluated and efficiently programmed in most 
standard software packages (e.g. MAPLE, MATHEMAT¬ 
IC A). We further demonstrate that previous results in (H, 
0 , a for the Rayleigh and Nakagami-m RF fading 
channels can be considered as special cases of our results. 

• We derive the generalized analytical expressions of the 
end-to-end outage probability and average bit error rate 
(BER) for the mixed RF/FSO dual-hop system. These 
exact formulations are more general, since they corre¬ 
spond to n-fi and 7 ^-/i fading in the RF link, and account 
for pointing errors, atmospheric turbulence and different 
modulation schemes in the FSO link. 

• The asymptotic closed-form outage probability in the 
high signal-to-noise ratio (SNR) regime is assessed, 
which can provide useful insights into the impact of 
system and fading parameters on the overall system 
performance. From an engineering perspective, we partic¬ 
ularly investigate the effect of pointing errors in the FSO 
link on the BER performance, over different modulation 
schemes and/or various channel fading conditions. 

• Although corresponding expressions for the case of n- 
/i/FF fading presented herein are given in terms of 
infinite series, we only need less than ten terms to get 
a satisfactory accuracy (e.g., smaller than 10 “^) for all 
considered cases. 

The remainder of the paper is structured as follows: In 
Section [II| the mixed RF/FSO system model is presented, and 
the statistical characteristics of generalized fading channels are 
described in Section Jill In Section JVl we derive novel and 
analytical expressions of the outage probability and average 
BER for the mixed RF/FSO dual-hop system. Numerical and 
Monte-Carlo simulation results are provided in Section [Vl and 
Section ED concludes with a summary of the main results. 

II. Mixed RF-FSO System Model 

As shown in Fig. [B we consider a dual-hop relay system 
with the source node S, the relay node R and the destination 


node D. The source node S cannot communicate directly with 
the destination node D because of the long distance or the 
block between them. The RF link between S and R experiences 
a /^-/i and 7 /-/i fading, since the ti-ji and r]-/! fading is more 
general and includes the Rayleigh, Rician and Nakagami-m 
fading as special cases ifTll . While the FSO link between 
R and D follows a FF distribution, which is widely used 
in FSO communications O. Similar to the related literature 
lEl-lIIl, we consider the high-energy FSO system with the 
commonly used on-off keying modulation. In this case, the 
FSO system performance is limited by background radiation 
and thermal noise, which can be modeled as the additive 
white Gaussian noise (AWGN) with the variance Nq/2 as a 
good approximation of the Poisson photon-counting detection 
model 1^ . Furthermore, we consider the typical amplify-and- 
forward (AF) relaying scheme in the relay node R, which 
simply amplifies its received signal by a fixed gain G and 
then forward it to the destination node D. Following the 
similar analysis in ED, the end-to-end received instantaneous 
electrical SNR 7 eq at the destination node D is given by 


7eq = 


7i72 
C + 72 ’ 


( 1 ) 


where c = G‘^/No, 71 denotes the SNR of the RF hop (i.e. 
S-R link) and 72 represents the electrical SNR of the FSO hop 
(i.e. R-D link). 


III. Statistical Characteristics 

In this section, we first deduce the novel analytical CDF 
expression for the mixed RF/FSO dual-hop system. Then, the 
PDF expression is derived by differentiating the CDF expres¬ 
sion. These expressions are useful for following analysis. 


A. Cumulative Distribution Function 


The CDF of the end-to-end instantaneous received electrical 
SNR 7 eq can be expressed as 


■ 7eq 


POO 

(x) = Pr [7eq < x\ 72]/^2 {72) ^72 

Jo 


C" ( 


c + 7 2 , 
72 


^] fl 2 (72)6^72. 


( 2 ) 


Note that ^ involves the PDF of the SNR at the relay node R, 
and the CDF of the electrical SNR at the destination node D. 
The PDF of FF fading channels with pointing errors is given 

by m, lEa 


/rr (72) 


tr (a)r(6)72 


^3,0 

^1,3 



e^ + 1 

( 3 ) 


where G(') is the Meijers’s G-function (241 Eq. (9.301)], r(-) 
is the gamma function (24l Eq. (8.310)], t accounts for the 
widely used detection technique type (i.e. t = 1 represents 
heterodyne detection, and t = 2 is intensity modulation with 
direct detection (IM/DD)), d = ^^/(^^ + 1) with ^ being 
the ratio between the equivalent beam radius and the pointing 
errors displacement standard deviation at the destination node 
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D, and Kt denotes the average SNR as 0, ll22ll 

'E( 72 )= 72 , t = l 

i2ahe{e+‘^) . „ (4) 


A 

l^t= i 


t = 2 




b = 


exp 


exp 


0.49(t| 




(l +0.18^2 + 0.56(72^/®)^'^^ 

0.51cr| (^1 + 0.69(72^/®) 


\ 


(l + O.QcP + 0.62(^2(72^/®) 


5/6 


- 1 


- 1 


r(/i) 

fi—l fi—k—l 


n=l 


l\ 


6:^2 ^ exp {nil) 


7i 


^ Ifi—i I 2/i 



n{l^ n)^ 


7i 


t (a+l) (6+l)(C2 + l)2’ 

Note that larger effect of pointing errors means smaller values 
of Moreover, the two effective numbers a and b are related 
to the atmospheric conditions and can be expressed as ini 

-1 


where /^(•) is the modified Bessel function of the first kind 
with order v. After using the infinite series representation of 
Iv{') EH Eq. (8.445)] and simple algebraic operation, the 
CDF of the n-fi distribution can be written as 

= (9) 


where A = 


A /x(l + Av) 


and r(',') is the incomplete gamma 


(5) 


, ( 6 ) 


function EH Eq. (8.350.2)]. 

For the case of ? 7 -/i/rr fading, we substitute ([3]) and (|7]) 
into El) and derive 

^2 




w here is the Rytov variance, d = 

kD‘^lAL with L being the distance between transmitter and 
receiver, k = 2 tt/\ denotes the optical wavenumber with A 
being the operational wavelength, D is the aperture diameter 
of the receiver, and is the altitude-dependent index of 
the refractive structure parameter determining the turbulence 
strength. Furthermore, we assume remains constant for 
a relatively long transmit bits interval, and it varies from 
1 X to 3 X for weak to strong 

turbulence cases EQl. 

Then we introduce the statistical characteristics of the gen¬ 
eralized fading for the RF hop. There are two formats for the 
r]-li distribution, namely. Format 1 and Format 2. According 
to imi, Format 1 can be converted into Format 2 through a 
bilinear transformation. Without loss of generality, we focus on 
Format 1 in the following. Due to its mathematical briefness 
and significant features ESI, the common case of integer 
values of is considered. The CDF of the 7 ^-/i distribution 
for integer values of /i is given by ESI Eq. (3)] 


tr (/i) r (a) r ( 6 ) \h 

2 fi—l fi—k—1 

xEE E 


^nQn,fc7Xxp(-74„7) 


n=lfe=0 (=0 


7 

X gI 


3,0 

3 


IH-exp 

72 


dab 


ll 

^7 \ 1 



+1 

V«t/ 

SY,a,b 


72 / 72 
^72. 


( 10 ) 


Utilizing the binomial expansion 1^ Eq. (1.111)] and the in¬ 
tegral identity as derived in (IA.4b . the corresponding analytical 
CDF of 7 eq can be derived as 

t2^a-\-b—2 f h\^ 

( 7 ) = 1_ — _( — 1 

r(/i)r(a)r(6)(27r)*-7^f/ 


2 ii—l^i—k—l I 

xEE E E 

n=l k=0 1=0 j=0 


^3t+l,0 


,k{Anl)^ ■’exp(-A„7) 


{dab)^AnC 

UJ 


T 


where uj 

|i! 

f t ’ 


a g+t-l h h+t-1 


r r+1 r+t 

^ t ^ t 




} and 


( 11 ) 


t ’ t ’ t ’ t ’ t 


!}■ 


xE E -^ 7 ^ 71 'exp (-.4„7 i) , (7) 

k=0 1=0 


where h = (2 + 77 ^ + 77 )/4, H = (77 ^— 77 ) 74 , 

Ai 4 2ii{h-H)lFi, A 2 ^ 2ii{hFH)/Fi, ai,fe ^ 

2 t^+kk\{h-H)F-^^ ^ 2 ,/c 2 i^+^k\{H+h)F-^ ’ denotes 

the average SNR of the RF hop. 

Moreover, the n-ii SNR PDF is given by Eq. (2)] 


( 7 .) = - exp f- " 


Eor the case of n-ii/TT fading, we utilize the finite series 
representation of r(-,-) EH Eq. (8.352.4)] and the similar 
method aforementioned. Then the corresponding analytical 
CDE of 7 eq can be derived as 


{kuY 


^2^a+6- 


F'^-^ ( 7 ) = - 

Fo (27r)*”7 (a) r ( 6 ) e^T'+«A‘ 

CO 74 +i 1 I / . \l —jf / \i 

^ (A 7 ) {nil) ^ 3 t+l,Q 

i\j\n_j\\ ^t,3t+i 
i=0 1=0 j=0 ^ "'z 


{dab)^Ac 

UJ 


T 


( 12 ) 


( 8 ) 


Although the CDE expression of n-ii/TT fading channels is 
given in the form of infinite series, we can truncated (fTH) 
suitably (no more than ten terms) so as to achieve a satisfactory 
accuracy (e.g., smaller than 10“^) as seen in Eig. O 


^Note that there is a typo for in (25] Eq. (3)], where {H - 
the denominator of ai^k should be {h — . 


in It is worthy to point out that for 77 -> 0, /i: ^ 0 and 11 = m 
(or 77 ^ 1 and /i = 77 i/ 2 ), ([TT]) and (fTH) reduce to the CDE 
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of the mixed Nakagami-m/Pr dual-hop relaying systems as 

—2 ^2 


^Nakagami ^ ^ 


Teq 

m—1 I 


(27r)‘-'r(a)r(6) 


X 


oi r/_oV I 




/o3t-|-l ,0 


m 

exp I ——7 
7i 

{dahf c m 

71 ^ 


CJ 

r 


which is in agreement with ||71 Eq. ( 8 )]. Furthermore, for 
dual-hop relaying systems over mixed Rayleigh/EE fading 
channels, e.g. 7 ^ 0 , and /i = 1 , our result presented 

in ([TT]) and ([T^ can be shown to agree with O Eq. ( 2 )]. 

Probability Density Function 

The PDF expression can be obtained by differentiating 
([TT1) and ([T^ with respect to 7 . With the help of (261 Eq* 
(07.34.20.0001.01)] and (24l Eq. (9.31.5)], we derive the 
analytical PDF expressions for two cases as 

p-U 

r(/i)r(a)r(6)(2^)'-^ V^y 

2 fi-l I r A \l 

n=l /c=0 Z=0 j=0 


JHi-jV- 


A„, - 


I-j 


^3t+l,0 


{dahf AnC 

UJ 

K.tC* ^ 



- ^3t+l,l 
■'-^t+l,3t+2 


{dahf AnC 


and 


•^7eq 


( 7 ) = 


^2 ^(2-|-6—2 


0, cc 

1 


(27r)"“"r(a)r(6) 

00 1 I / . \Z—i / \i 

(^ 7 ) (^/^) 


^ ^ ^ ^ 

Z^ o\n\ 


ilj! (/-?■)! 
i=o z=o i=o ^ 


^3t+l,0 


Z- 


l-j 


{dab)^Ac 

UJ 

KtP* 

rj 


_ ^^3t+l,l 
^W+l,3t+2 


{dabfAc 


0,u; 
r, 1 


probability of dual-hop AF relaying can be obtained by setting 
7 = 7th in (HU) as Po = Pr ( 7 eq < 7 *) = P^,, ( 7 *). 

In order to provide more insights into the impact of system 
and fading parameters on the system performance, we now 
focus on the high-SNR regime. By taking Kt ^ I in ([TT]) and 
(O, and using (24l Eq. (9.303)], the asymptotic high-SNR 
outage probability can be derived as 




^2^a-|-6 


-2 


r(/x)r(a)r(6 ) (27r) 

2 /X—1 /X—/c—1 I 3t+l 


( p, 


<^n,/c (^n7th) 


1-3 


EE E EE 'll il — 7)^6 -^nTth 

n=l/c=0 1=0 j=0 x=l 


/ (dab) AnC 


UU^H-n) ’ 


(15) 


and 


pK-ll 


_ _ 

Z—^ 

i=0 

00 Ai-|-x—1 I 


^2^a+b 


-2 


E 


(27r)*”^r (a) r ( 6 ) e^'r.h+«^ 

(A7,h)'“7«M)* 


i=o /=o i=o 
/ (dab)* Ac 

X I -^;^7th 






3t+l 

p=l,p^i 




■i) 


Ktfit '7 ]yp=i^{^p-Ti) 


(16) 


(13) 


(14) 


where 7 denotes the i-th term of r, and ujp accounts for the 
p-th term of cc, respectively. Note that (fTSl) and ([16]) are in the 
form of simple elementary functions. It is clear to see from 
([T5]) and (O that the outage probability becomes smaller as 
the average electrical SNR Ht of the FSO link increases, which 
agrees with the corresponding findings presented in iffl-llTl. 

B. Bit Error Rate 

Note that the binary modulation schemes are widely used 
in experimental and practical FSO communication systems 
(28l, ||29]| . For an extensive list of different binary modulation 
schemes, a unified BER expression is given as (301 Eq. (12)] 

qP foo 

-Pb = 77FT77T / exp (-^ 7 ) 7 P"^P ( 7 )(i 7 , (17) 


Note that ([T^ and (El) can be reduced to m Eq. (11)] and 
(S] Eq. (4)] for the cases of Nakagami-TTT'/FF and Rayleigh/FF 
fading, where the IM/DD scheme is considered, respectively. 

IV. Performance Analysis 

We derive novel and analytical expressions of the end-to-end 
outage probability and average BER for the mixed RF/FSO 
dual-hop system in this section. Moreover, the asymptotic out¬ 
age probability in the high-SNR regime and BER performance 
without pointing errors are also assessed. 

A. Outage Probability 

As an significant metric for the performance of a wireless 
communication system, the outage probability is defined as 
the probability that the instantaneous SNR 7 eq falls below a 
predetermined threshold 7 th (271. Mathematically, the outage 


2r (p) Jo 

where the parameters p and q denote different modulation 
schemes. For example, by setting p = 1 and q = 0.5, the 
BER of non-coherent binary frequency shift keying (NBFSK) 
is obtained, while for coherent binary frequency shift keying 
(CBFSK), we should set p = 0.5 and q = 0.5 (til- 
By substituting ([TT]) and ([T2b into ([TT]) . and using (24l Eq. 
(7.813.1)], we can obtain the analytical BER expressions as 

1 p^a+h-2^p ■ 

2 ~ 2r7j 

2 p-1 p-k-1 I \l-j / , A \j-p-l 




xEE E E 

n=l k=0 1=0 j=0 


<^n,/c(^n) \q F An) 

JHi-JV- 


X G 


-3t+l,l 

t+l,3t+l 


{dabY AnC 


KtC* {q + An) 

T 


( 18 ) 














































5 


and 


pK-tl _ 1 ^ 1^^) 


^2^a+b—2^p^—KH 


2^i!e'=M 2r(Ai)r(a)r(6)r(p)(27r) 


t-1 


i=0 Z=0 j=0 


i\j\{i-j)\ 


X G: 


3t+l,l 

t+l,3t+l 


(dab) Ac 


{q + 


l+j-p-l,uj 

T 


(19) 

Note that for the special case of mixed Rayleigh/Pr fading 
channels as investigated in (Si, (fTSl) and (O can reduce to (H 
Eq. (14)] where the IM/DD scheme is considered. 

Then, we consider the non-pointing errors case and recall 
that d = ^^/(^^ + 1). By taking ^ ^ oo (larger ^ means 
smaller pointing errors) and using the definition of the Meijer’s 
G-function ||24l Eq. (9.301)], the BER expressions can be 
converged to 

^a+b-2^p ^ 


1 




2r (Ai)r(a)r(6)r(p) (27r) 


t-1 


H 


2 fi—lfi—k—1 I 

xEE E E 

n=l k=0 1=0 j=0 


'-^t-l,2t+l 


(^n,k{Afi^ ip H“ Aji) 

MET 


j-p-l 


{abf AnC 

^+j-p-l 

{q +An) 

t' 

( 


( 20 ) 


and 

rthi-p 


9 2^ il 


^2^a-\-b-2qP^-Kp 


b.«^oo 2 ^ 2r (/x) r (a) T (b) T (p) (27r) 


t-1 


t=0 Z=0 ji'=0 


^ '-^t-l,2t+l 


z!j!(;-j)! 


{ah)* AnC 

1 + i - p - ( 

KtC* {q +An) 

t' 


( 


( 21 ) 


where r' = {j, k±LA^jy Once again, it is worthy 

to point out that (l20l) and (1211) reduce to |[5] Eq. (15)] when t = 
2 (i.e. IM/DD scheme). Note that the average BER expressions 
of non-binary modulation schemes can be also obtained by 
using the similar way presented here. 


V. Numerical Results 

In this section, the analytical results of mixed RE/ESO 
systems over different system configuration scenarios and/or 
various channel fading conditions are presented and compared 
with Monte-Carlo simulations. Without loss of generality, 
we assume the heterodyne detection t = 1, the threshold 
7th = OdB, the link distance L = 4000m, the receiver aperture 
diameter D = 0.01m, and a wavelength of A = 1550nm 
throughout this section 1^ . Eor the mixed RE/ESO AE relay 
scheme with the fixed gain, we set c = 1 at the relay node. 
In Monte-Carlo simulations, 10^ r]-p and EE random samples 
are generated. 



Fig. 2. Simulated, analytical and high-SNR approximated outage probability 
of mixed RE/ESO dual-hop relay systems over r]-/j,/rr fading channels versus 
the average electrical SNR of the FSO hop 72 (C^ = 1 x 10“ ^^m“^/^, 
D = 0.01m, L = 4000m, A = 1550nm, 71 = lOdB and ^ = 1.1). 



Fig. 3. Simulated, analytical and high-SNR approximated outage probability 
of mixed RE/ESO dual-hop relay systems over K-j^/FF fading channels versus 
the average electrical SNR of the FSO hop 72 (C^ = 1 x 10“ ^^m“^/^, 
D = 0.01m, L = 4000m, A = 1550nm, 71 = lOdB and ^ = 1.1). 


Eigure [2| depicts the simulated, analytical (fTTI) and high- 
SNR approximated (fTSl) outage probability against the average 
electrical SNR of the ESO hop when different values of the 
fading parameters /i and r] are considered in the RE hop. The 
average SNR of the RE hop is set as = lOdB, the turbulence 
strength is weak with G^ = 1 x 10 “and the pointing 
errors are large with ^ = 1.1. We can find that there is a 
good match between the analytical and simulated results, and 
the high-SNR approximation is quite tight in the moderate- 
and high-SNR regimes, which validates the accuracy of the 
proposed expressions. By varying one parameter while keeping 
the other one fixed, Eig. [2] reveals that increasing the values 
of fi and/or r] helps to overcome the fading effects. This 
is due to the fact that with more multipath clusters (larger 
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Fig. 4. Simulated and analytical outage probability of mixed RF/FSO dual¬ 
hop relay systems over weak to strong turbulence fading channels versus the 
average electrical SNR of the RF hop = 3, r] = 0.5, D = 0.01m, 

L = 4000m, A = 1550nm, 7-2 = lOdB and ^ = 1.1). 



Fig. 5. Simulated and analytical BER of mixed RF/FSO dual-hop relay 
systems over r]-fi and moderate turbulence fading channels versus the average 
electrical SNR of the RF hop 'fi (iJ. = 3, r] = 0.5, D = 0.01m, L = 4000m, 
A = 1550nm, ^2 = lOdB and = 9 x lO-^^m-^/S). 


/i) and/or the less difference between the power of in-phase 
and quadrature components (larger r]), the power of received 
signals will be enhanced. We recall that consistent conclusions 
were also drawn in EH, m. 

It is clear to see from Fig. |3] that the simulated, analytical 
(Ell and high-SNR approximated (O outage probability 
results match well. This verifies the accuracy of the CDF 
expression (Ell derived in the analysis. With the selected 
parameters as shown in Fig. [3l we find that the convergent 
series of ([T^ is truncated with only ten terms to be sufficient 
to get numerically accurate results. In the tz-fi model, the 
parameter k. indicates the power of the dominant component 
of a signal and the parameter fi represents the number of 
multipath clusters, respectively. By varying one parameter 


while keeping the other parameter fixed. Fig. [3] reveals that 
increasing the values of or /i both helps overcome the effects 
of fading. Moreover, the effect of fi on the outage probability is 
more pronounced than that of /^. For example, the gap between 
/i: = 3, /i = 1 and /i: = 3, /i = 2 curves is larger than that 
between /i: = 3, /i = 1 and k, = 3, jj. = 2 curves. 

In Fig. m the outage probability of mixed RF/FSO dual-hop 
relay systems under rj-fi/TT fading conditions with 77 = 0.5 
and /i = 3, and strong pointing errors with ^ = 1.1 
is illustrated for various values of the strength of the at¬ 
mospheric turbulence C^. The infiuence of the atmospheric 
turbulence on the outage performance becomes more obvious 
from weak to strong turbulence channels. For example, at 
7 i = 20dB, the difference of the outage probability between 
= 9 X 10“^^m“^/^ and = 1 x 10“^^m“^/^ is 

much larger than that between = 3 x 10 “^^m“^/^ and 

= 9 X 10“^^m“^/^. Therefore, Fig. [4] indicates that the 
mixed RF/FSO dual-hop relay system losses much outage per¬ 
formance in adverse weather or strong turbulence conditions 
(i.e. large values of C^). 

The effect of pointing errors on the BER performance of 
mixed RF/FSO dual-hop relay systems with different modu¬ 
lation schemes is studied in Fig. [Jl The parameters of the RF 
link are r] = 0.5 and /i = 3, and the turbulence strength is 
= 9 X 10“^^m“^/^. It is clear that the analytical BER 
expressions coincide with the simulation results. We can also 
find from Fig. [5] that the BER of CBFSK is smaller than 
that of NBFSK. This observation can be explained by the 
fact that the system achieves better BER performance when 
the channel state information is available at the receivers. 
More importantly, the pointing errors (note that larger ^ means 
smaller pointing errors) have an obvious effect on the BER 
performance, especially when ^ is small. This observation 
agrees with the corresponding finding in 0. 

VI. Conclusions 

In this paper, we investigated a dual-hop relaying system 
over mixed RE/ESO links, which can be modeled as 77 - 
/i/FF and /^-/i/FF distributions. Novel and exact expressions 
of the CDF and PDF of the end-to-end SNR were derived 
in terms of finite sums of Meijer’s G-functions. We also 
derive the new analytical expressions of the end-to-end outage 
probability and average BER for such system. These results 
are general, since they correspond to n-ii and r]-ii fading 
in the RE link, and account for pointing errors, atmospheric 
turbulence and different modulation schemes in the ESO link. 
Note that for the case of /^-/i/FF fading, the infinite series in 
corresponding results quickly and steadily converges, requiring 
only a few terms to obtain a desired accuracy. For example, 
to obtain an error smaller than 10“^, less than ten terms 
are required in all of considered cases. Capitalizing on the 
high-SNR asymptotic outage probability, some useful insights 
were revealed. Particularly, it is clear to find that the outage 
probability becomes smaller as the electrical SNR of the FSO 
link increases. In addition, we demonstrated that the pointing 
errors have an obvious effect on the BER performance. 
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Appendix 


We involve the integral in terms of power, exponential and 
Meijer’s G-functions as 

poo 

Jo 

(A.l) 


y,u/v 


(cLp) 

(b,) 


dx. 


where (a^) = {ai, • • • , a^}, {hq) = {6i, • • • , hq). Applying a 
change of variables, x~^ ^ x and with the help of ||24| Eq. 
(9.31.2)], (lA.ll) becomes 


pOO 

= / exp (-aa;)G”’™ ^ 

Jo ^ 9.P L l-(ap) 


1 - {b,) 


v°u 


a-1/2 -a 


(27r) 


(w-l)/2+c*(v-l) 


p-^vn,vm-\-u 


dx 

(A.2) 


UJ '^u'^ 

A{u,l-a),A{v,l- (bg)) 

^Uyv{p-q) 

A{v,l- (ap)) 


(A.3) 


where s = E (1 “ ®i) ~ E (1 ~ bj) + + 1, c* = m + 

7=1 7=1 

n- (p + g) /2, and A {u,v) = • • • , Note that 

from (lAAb to (IA.3K we have used the identity of 11331 Eq. 
(2.24.3.1)]. By using ||24l Eq. (9.31.2)] again, the analytical 
result of (lA.lb can be expressed as 


/ = 


r^jSyCX 1/2^ a 


X G 


vm+u,vn 

vp^vq\u 



l-A{v,l-{ap)) 

y^UyV{q-p)) 

l-A{u,l-a),l- A{v,l-{bq)) 


(A.4) 


References 

[1] J. A. Simpson, B. L. Hughes, and J. F. Muth, “Smart transmitters and 
receivers for underwater free-space optical communication,” IEEE J. Sel. 
Areas Commun., voL 30, no. 5, pp. 964-974, June 2012. 

[2] Z. Ghassemlooy, W. Popoola, and S. Rajbhandari, Optical wireless 
communications: System and channel modelling with Matlab®. CRC 
Press, 2012. 

[3] N. D. Chatzidiamantis, H. G. Sandalidis, G. K. Karagiannidis, and 
M. Matthaiou, “Inverse Gaussian modeling of turbulence-induced fading 
in free-space optical systems,” J. of Lightw. Technol, vol. 29, no. 10, 
pp. 1590-1596, May 2011. 

[4] E. Lee, J. Park, D. Han, and G. Yoon, “Performance analysis of the 
asymmetric dual-hop relay transmission with mixed RF/FSO links,” 
IEEE Phot. Technol. Lett., vol. 23, no. 21, pp. 1642-1644, Nov. 2011. 

[5] I. S. Ansari, F. Yilmaz, and M.-S. Alouini, “Impact of pointing errors 
on the performance of mixed RF/FSO dual-hop transmission systems,” 
IEEE Wireless Commun. Lett., vol. 2, no. 3, pp. 351-354, June 2013. 

[6] H. Samimi and M. Uysal, “End-to-end performance of mixed RF/FSO 
transmission systems,” lEEE/OSA J. Opt. Commun. Netw., vol. 5, no. 11, 
pp. 1139-1144, Nov. 2013. 

[7] E. Zedini, I. Ansari, and M.-S. Alouini, “Performance analysis of mixed 
Nakagami-m and Gamma-Gamma dual-hop FSO transmission systems,” 
IEEE Photon. J., vol. 7, no. 1, pp. 1-20, Feb. 2015. 

[8] E. Soleimani-Nasab and M. Uysal, “Generalized performance analysis 
of mixed RF/FSO systems,” in Proc. IEEE Opt. Wireless Commun. 
(IWOW), Funchal, Portugal, Sept. 2014, pp. 16-20. 

[9] M. A. Kashani, M. Uysal, and M. Kavehrad, “A novel statistical model 
for turbulence-induced fading in free-space optical systems,” in Proc. Int. 
Conf. Transparent Opt. Netw. (ICTON), Cartagena, Spain, Jun. 2013, pp. 
1-5. 

[10] J. D. Parsons, The mobile radio propagation channel, vol. 2. John Wiley 
New York, 2000. 


[ 11 ] 

[ 12 ] 

[13] 

[14] 

[15] 

[16] 

[17] 

[18] 

[19] 

[ 20 ] 

[ 21 ] 

[ 22 ] 

[23] 

[24] 

[25] 

[26] 

[27] 

[28] 

[29] 

[30] 

[31] 

[32] 

[33] 


M. D. Yacoub, “The n-p, distribution and the p-p distribution,” IEEE 
Antennas Propag. Mag., vol. 49, no. 1, pp. 68-81, Feb. 2007. 

L. C. Andrews and R. L. Phillips, Laser beam propagation through 
random media. SPIE press, 2005. 

N. Wang and J. Cheng, “Moment-based estimation for the shape 
parameters of the Gamma-Gamma atmospheric turbulence model,” Opt. 
Express, vol. 18, no. 12, pp. 12 824-12 831, June 2010. 

M. A. Al-Habash, R. L. Phillips, and L. C. Andrews, “Mathematical 
model for the irradiance probability density function of a laser beam 
propagating through turbulent media,” Opt. Eng., vol. 40, no. 8, pp. 
1554-1562, Aug. 2001. 

A. Jurado-Navas, J. M. Garrido-Balsells, J. F. Paris, M. Castillo- 
Vazquez, and A. Puerta-Notario, “Impact of pointing errors on the 
performance of generalized atmospheric optical channels,” Opt. Express, 
vol. 20, no. 11, pp. 12 550-12 562, 2012. 

A. Garcia-Zambrana, C. Castillo-Vazquez, B. Castillo-Vazquez, and 

R. Boluda-Ruiz, “Bit detect and forward relaying for FSO links us¬ 
ing equal gain combining over gamma-gamma atmospheric turbulence 
channels with pointing errors,” Opt. Express, vol. 20, no. 15, pp. 16 394- 
16409, 2012. 

E. Bayaki, R. Schober, and R. K. Mallik, “Performance analysis of 
MIMO free-space optical systems in gamma-gamma fading,” IEEE 
Trans. Commun., vol. 57, no. 11, pp. 3415-3424, Nov. 2009. 

N. D. Chatzidiamantis and G. K. Karagiannidis, “On the distribution of 
the sum of gamma-gamma variates and applications in RF and optical 
wireless communications,” IEEE Trans. Commun., vol. 59, no. 5, pp. 
1298-1308, May 2011. 

M. Aggarwal, P. Garg, and P. Puri, “Dual-hop optical wireless relaying 
over turbulence channels with pointing error impairments,” J. Lightwave 
Technol, vol. 32, no. 9, pp. 1821-1828, May 2014. 

S. M. Navidpour, M. Uysal, and M. Kavehrad, “BER performance 
of free-space optical transmission with spatial diversity,” IEEE Trans. 
Wireless Commun., vol. 6, no. 8, pp. 2813-2819, Aug. 2007. 

K. P. Peppas, A. N. Stassinakis, H. E. Nistazakis, and G. S. Tombras, 
“Capacity analysis of dual amplify-and-forward relayed free-space op¬ 
tical communication systems over turbulence channels with pointing 
errors,” IEEE/0SA J. Opt. Commun. Netw., vol. 5, no. 9, pp. 1032-1042, 
Sep. 2013. 

H. G. Sandalidis, T. A. Tsiftsis, and G. K. Karagiannidis, “Optical 
wireless communications with heterodyne detection over turbulence 
channels with pointing errors,” J. Lightw. Technol, vol. 27, no. 20, pp. 
4440-4445, Oct. 2009. 

W. Gappmair, “Further results on the capacity of free-space optical 
channels in turbulent atmosphere,” lET Commun., vol. 5, no. 9, pp. 
1262-1267, May 2011. 

I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and 
Products, 7th ed. San Diego, CA: Academic Press, 2007. 

K. P. Peppas, “Dual-hop relaying communications with cochannel inter¬ 
ference over p-p fading channels,” IEEE Trans. Veh. Technol, vol. 62, 
no. 8, pp. 4110-4116, Oct. 2013. 

Wolfram, “The Wolfram functions site,” Available: 
http://functions.wolfram.com. 

M. K. Simon and M.-S. Alouini, Digital communication over fading 
channels. Jone Wiley & Sons, 2005. 

W. O. Popoola and Z. Ghassemlooy, “BPSK subcarrier intensity modu¬ 
lated free-space optical communications in atmospheric turbulence,” J. 
Lightwave Technol, vol. 27, no. 8, pp. 967-973, april 2009. 

M. Karimi and M. Nasiri-kenari, “BER analysis of cooperative systems 
in free-space optical networks,” J. Lightwave Technol, vol. 27, no. 24, 
pp. 5639-5647, Dec 2009. 

I. S. Ansari, S. Al-Ahmadi, F. Yilmaz, M.-S. Alouini, and 

H. Yanikomeroglu, “A new formula for the BER of binary modulations 
with dual-branch selection over generalized-A:,” IEEE Trans. Commun., 
vol. 59, no. 10, pp. 2654-2658, Oct. 2011. 

I. Trigui, A. Laourine, S. Affes, and A. Stephenne, “Performance 
analysis of mobile radio systems over composite fading/shadowing 
channels with co-located interference,” IEEE Trans. Wireless Commun., 
vol. 8, no. 7, pp. 3448-3453, July 2009. 

J. Zhang, M. Matthaiou, Z. Tan, and H. Wang, “Performance analysis 
of digital communication systems over composite ? 7 -/i/gamma fading 
channels,” IEEE Trans. Veh. Technol, vol. 61, no. 7, pp. 3114-3124, 
Sep. 2012. 

A. P. Prudnikov, I. U. A. Brychkov, and O. I. Marichev, Integrals and 
Series, Volume 3: More Special Eunctions. Gordon and Breach, 1990. 



















